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MULTI-COMPARTMENT SEPARATION IN MAGNETIC RESONANCE USING 
TRANSIENT STEADY-STATE FREE PRECESSION IMAGING 



GOVERNMENT RIGHTS 

[0001] The U.S. government has rights in the disclosed invention pursuant to NIH grants 
to Stanford University including NIH-HL39297, NIH-HL56394, NEH-AR46904, and NIH- 
CA50948. 



BACKGROUND OF THE INVENTION 

[0002] This invention relates generally to magnetic resonance imaging (MRI) of tissue 
having a plurality of tissue species with different relaxation times, and more particularly the 
invention relates to separation of the species from analysis of MRI relaxation curves using . 
transient steady-state free precession (balanced SSFP) imaging. 

[0003] Magnetic resonance imaging (MRI) is a widely used medical imaging modality 
that provides excellent soft-tissue contrast with arbitrary scan-volume orientations. Unlike 
X-ray computed-tomography or ultrasound, whose contrast is based only on the transmission 
or reflection properties of tissue, MRI generates contrast from a variety of physical properties 
of tissues including relaxation, chemical-shift, diffusion and proton density. 
[0004] Clinically, one of the most useful MRI contrast mechanisms is T 2 -contrast, which 
arises from differences in the spin-spin relaxation time (T 2 ). T 2 -contrast has numerous uses 
including distinguishing malignant and benign tumors, imaging spinal abnormalities, vascular 
imaging and for diagnosing meniscal tears in the knee. Quantitative measurement of the T 2 
spectrum, or T 2 -relaxometry, is also useful for a variety of clinical applications including 
cartilage imaging and imaging of multiple sclerosis (MS). 

[0005] T 2 -relaxometry is discussed by Harrison et al. in "Magnetization Transfer and T 2 - 
Relaxation Components in Tissue," Magnetic Resonance in Medicine 1995; 33: 490-496. 
Many tissues are known to exhibit multi-component T 2 -relaxation that suggests some 
compartmental segregation of mobile protons on a T 2 time scale. Magnetization transfer 
(MT) is another relaxation mechanism that can be used to produce tissue contrast in MR 
imaging. The MT process depends strongly on water-macro molecular interactions. To 
investigate the relationship between multi-component T 2 -relaxation and the MT process, 
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multi-echo T2 measurements have been combined with MT measurements. For example, in 
muscle, short-T 2 components show greater MT than long-T 2 components while for white 
matter, MT measurements were identical for two major T 2 components, apparently because 
of exchange between the T 2 compartments on a time scale characteristic of the MT 
experiment. 

[0006] Whittal and Mackay, "Quantitative Interpretation of NMR Relaxation Data," 
Journal of Magnetic Resonance 84, 134-152 (1989), discloses a least-squares and linear 
programming algorithms for the interpretation of NMR relaxation data in terms of a spectrum 
of relaxation times. 

[0007] In the case of MS, the exact T 2 -spectrum is not as important as the fraction of the 
spectrum that appears in each peak. Measurement of the relative peak sizes is a good 
indicator of the myelination of white matter in the brain, and an indicator of MS. Currently, 
this measurement is performed using a multi-echo Carr-Purcell Meiboom-Gill (CPMG) 
sequence, see Whittal and Mackay, supra. In the CPMG sequence, the minimum echo 
spacing is on the order of 10-15 ms to allow for crusher gradients and 180° refocusing RF 
pulses. 

[0008] The present invention presents a new method for quantitatively separating tissue 
based on relaxation time differences with shorter repetition time and improved temporal 
resolution. 

SUMMARY OF THE INVENTION 

[0009] In accordance with the invention, the exponential transient signal decay of a 
steady-state free precession (balanced SSFP) imaging sequence is utilized in obtaining multi- 
compartment separation of different tissue species, similar to the standard CPMG method of 
T 2 -relaxometry. However, the short repetition time of balanced SSFP (around 3-5 ms) allows 
time samples on the transient curve to be 3-4 times as closely spaced as the CPMG sequences 
allow. Further, the exponential decay rate of balanced SSFP is slower than T 2 thereby 
allowing more time samples before the signal from short-T 2 species has decayed. 
[0010] The invention and objects and features thereof will be more readily apparent from 
the following detailed description and appended claims when taken with the drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Fig. 1 illustrates a steady-state free precession imaging sequence preceded by a 
series of catalyzation flip pulses where acquisition begins immediately following the 
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catalyzation period and the phase of each RF pulse alternates by 180° on each sequence 
repetition. 

[0012] Fig. 2 illustrates brain images during a transient signal period as the steady-state 
magnetization evolves where the image signal is a combination of exponentially-decaying 
signals from cerebral-spinal fluid (dash-dot line) extracellular water in white matter (solid 
line) and water in myelin (dashed line). 

[0013] Fig. 3 illustrates r- spectra obtained from a series of transient images shown in Fig. 
2 where for a region of interest, the total signal of each image is fitted using a non-negative 
curve fit to determine the t spectra and with the separate peaks in the t spectra corresponding 
to myelin water and extracellular water in white matter. 

DETAILED DESCRIPTION OF ILLUSTRATED EMBODIMENTS 
[0014] The use of balanced SSFP imaging is advantageous compared with CPMG for short- 
T 2 species for two reasons. First, the short repetition time of balanced SSFP, around 3-5 ms 
allows the time samples on the exponential curve that are 3-4 times as closely spaced as the 
CPMG sequences allow. Second, the exponential decay rate is slower than T 2 , allowing more 
time samples before the signal from short-T2 species has decayed. 

[0015] Standard balanced SSFP imaging consists of RF excitation pulses spaced apart by 
the sequence repetition time, TR. All imaging gradients are rewound, and the low spatial 
frequency information is acquired at an echo time (TE) midway between RF excitation 
pulses. A steady-state signal arises after many repetitions. The steady-state signal provides a 
very efficient imaging method as well as useful contrast. 

[0016] In balanced SSFP imaging, magnetization gradually approaches a steady state 
when the signal is the same from one repetition to the next. During the time in which the 
steady state evolves, the transient signal can be well-predicted if the starting magnetization is 
known. Specifically, if the starting magnetization is roughly a scalar multiple of the steady- 
state magnetization, then the transient signal decays along a smooth exponential path to the 
steady state. The magnetization can be manipulated roughly to its steady-state direction 
using some smoothly increasing series of flip angles toward the flip angle of the balanced 
SSFP sequence. Data acquisition begins following this "catalyzation" period, as shown in 
Fig. 1. 

[0017] Images are acquired by repeating the sequence shown in Fig. 1 once for every 
phase encode line that needs to be acquired to make an image. This is similar to the Tj- 
measurement method presented by Scheffler and Hennig, "Ti Quantification with Inversion 
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Recovery TrueFISP," Magn Reson Med 2001; 45: 720-723, except that no inversion pulse is 
used, a larger flip angle is used, and the mixing of Ti and T2 is not a problem. 
[0018] A sufficiently high number of sequence repetitions, such as 256, is used to sample 
the exponential decay long enough that the signal from voxels of interest reaches a steady 
state. Following this, all RF and gradient pulses are turned off for a recovery period of about 
one second, and then the sequence is repeated for the next phase encode line. In preliminary 
trials, a TR of 5 ms and a flip angle of 60° were used. 

[0019] The result of this acquisition is a series of complete images of the transient decay, 
spaced apart by the repetition time of the balanced SSFP sequence. A representative example 
is shown in Fig. 2. The figure shows the exponential responses that should arise from the 
short-T2 myelin-water compartment, the extracellular water in white matter and the cerebral- 
spinal fluid (CSF). The observed signal decay in a voxel is a linear combination of decaying 
exponentials such as those shown. 

[0020] A small region-of-interest (ROI) can be selected in the image, as shown in Fig. 3. 
The sum of all pixels in this ROI is plotted as a function of time in the bottom-left of Fig. 3. 
The curve was fitted to a model: 

A*(0 = £ C f* * +M SS (1) 
1=1 

where M{t) is the observed signal intensity as a function of time, N is the number of 7 
points used in the fit, Q is the relative amount of material with an exponential time constant 
7. and M ss is the steady-state signal in the voxel. 

[0021] The coefficients C, were determined using a non-negatively constrained LI -norm 
method (similar to Xu et al., "Homogeneous Magnet Design Using Linear Programming," 
IEEE Trans Med Imaging 2000; 36: 476-483) for 7 values spaced apart by 30 ms. The 
resulting "r-spectra" are shown in the bottom-right plot in Fig. 3. The two peaks in the t- 
spectra correspond to the myelin-water and extracellular water as previously reported in 
Mackay et al., "In vivo Visualization of Myelin Water in Brain by Magnetic Resonance," 
Magnetic Resonance in Medicine 1994; 31 :673-677. 

[0022] T 2 -relaxometry has been performed by numerous investigators for many different 
purposes. The most common pulse sequences used for quantitative T 2 measurements are 
single-echo and multiple-echo (CPMG) spin-echo sequences (see Harrison et al. and Whittal 
and Mackay, supra). Single-echo methods are more accurate, but are very slow. CPMG 
methods have a minimum echo spacing of 10-15 ms depending on gradient performance and 
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resolution. For short-T 2 species, this means only a few samples are obtained on the decay 
curve before the signal is diminished to the level of background noise. 
[0023] The method in accordance with the invention has two primary advantages over 
CPMG techniques. First, the time samples can be spaced very closely. For image resolution 
of 1 mm, on current MR systems, the spacing can be as low as 3-5 ms, depending on the 
imaging trajectory used. For lower resolution images, the spacing could be reduced further. 
[0024] The second advantage of the current method is that the decay constant r is a 
combination of Ti and T 2 . When the goal of the imaging is to separate different 
compartments, such as myelin water and extracellular water in white matter, the mixing of Ti 
and T 2 still allows separation. However, the fact that the decay constant r is greater than T 2 
means that the decay of short-T 2 , species such as myelin water, is slowed. This generally 
reduces the error in estimating the amount of myelin water. 

[0025] The combination of these two advantages means that about 10 times as many 
samples can be obtained on the decay curve for short-T 2 species such as myelin water. This 
can improve the reliability of the curve-fit even though the signal level is lower than in a 
CPMG sequence. Further, there are 2-4 more samples along the curve. Additionally, flip 
angle can be used to manipulate the decay rate. 

[0026] The sampling method for multi-compartment analysis has been described using a 
simple Cartesian imaging readout and a standard balanced SSFP sequence with the 
catalyzation scheme described in Deimling and Heid, "Magnetization Prepared True FISP 
Imaging," In: Proceedings of the 2 nd Annual Meeting of SMR, San Francisco, 1994. p. 495. 
However, for multi-compartment analysis, the same phase-encoded frame is acquired for the 
entire transient decay, rather than acquiring multiple frames on one decay. The next phase- 
encoded frame is acquired over a subsequent transient decay. There are several variations of 
this method that may be useful. 

[0027] First, the Cartesian imaging trajectory is slow. A spiral or echo-planar imaging 
trajectory can reduce the total imaging time by a factor of about 3 for balanced SSFP imaging. 
Alternatively, to reduce imaging time, some interleaving could be used by trading off some of 
the high sampling along the exponential trajectories that this method provides. 
[0028] The variation of flip angle across the slice profile may cause problems. However, 
the effect of slice profile can be significantly reduced by using a 3D imaging trajectory. This 
is particularly useful in cases where signal averaging would otherwise be necessary to 
improve SNR. 

[0029] Another variation of this sequence is to precede the magnetization preparation 
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with an inversion pulse as described by Scheffler and Hennig, "Ti Quantification with 
Inversion Recovery TrueFISP," Magn Reson Med 2001 ; 45: 720-723. This effectively 
changes the values of the constants CJ, increasing them in magnitude and negating them. 
The constrained fit would be changed accordingly to force non-positive C_i. 
[0030] The invention realizes a higher achievable sampling density for relaxometry, as 
well as the ability to vary the decay times of different species for measurement. Particularly 
for short-T2 materials, both of these features should improve the accuracy and the robustness 
of curve-fitting the exponential decays. The invention utilizes a standard balanced SSFP 
(FIESTA, TrueFISP, or balanced-FFE) imaging sequence and a preparation sequence, usually 
consisting of RF tips identical to those used in the balanced SSFP sequence, but of varying 
amplitude and/or phase and spaced apart by TR or TR/2. A sufficiently short repetition time 
(TR) is used such that the time constant r does not vary much due to off resonance effects. 
Additionally, a curve-fitting algorithm such as non-negative least-squares (see Whittal and 
Mackay, supra) is used to fit multi-exponential decay. 

[0031] There has been described a new method of quantitatively separating tissue based 
on relaxation time differences. The method uses the transient signal decay in balanced SSFP 
imaging to provide an alternative to the standard CPMG methods of T 2 -relaxometry. The 
balanced SSFP technique allows 3-4 times the temporal resolution of CPMG, and also slows 
the short-T2 decay so that it can be more accurately measured. 

[0032] While the invention has been described with reference to specific embodiments, 
the description is illustrative of the invention and is not to be construed as limiting the 
invention. For example, a simple data frame can be acquired repeatedly over the entire decay 
curve of the balanced SSFP imaging sequence or multiple data frames can be acquired 
repeatedly over the decay curve. Further, an inversion pulse can be applied in the preparation 
stage whereby magnetization starts at a negative value. Alternatively, magnetization can be 
saturated in the preparation stage with magnetization starting at a zero value, or 
magnetization can be inverted in the steady state with the magnetization signal starting at a 
negative value. Various modifications and applications may occur to those skilled in the art 
without departing from the true spirit and scope of the invention as defined by the appended 
claims. 
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